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The goal of the present paper is the investigation of the evolution of anisotropic regular structures 
and turbulence at large Reynolds number in the multi-dimensional Burgers equation. We show that 
we have local isotropization of the velocity and potential fields at small scale inside cellular zones. 
For periodic waves, we have simple decay inside of a frozen structure. The global structure at 
large times is determined by the initial correlations, and for short range correlated field, we have 
isotropization of turbulence. The other limit we consider is the final behavior of the field, when 
the processes of nonlinear and harmonic interactions are frozen, and the evolution of the field is 
determined only by the linear dissipation. 
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I. INTRODUCTION 

The well known Burgers equation describes a variety 
of nonlinear wave phenomena arising in the theory of 
wave propagatioUjacoustics, plasma physics and so on 
(see, e.g., [Illl, [1,11, [1,0]) ■ This equation was originally 
introduced by J. M. Burgers as a model of hydrodynam- 
ical turbulence [3, [1] ■ It shares a number of properties 
with the Navier-Stokes equation : the same type of non- 
linearity, of invariance groups and of energy-dissipation 
relation, the existence of a multidimensional version, etc 
0. However, Burgers equation is known to be integrable 
and therefore lacks the property of sensitive dependence 
on the initial conditions. Nevertheless, the differences 
between the Burgers and Navier-Stokes equations are as 
interesting as the similarities ^Td\ and this is also true for 
the multi-dimensional Burgers equation : 



'dt 



+ (v • V) V = i^V^v . 



(1) 



With external random forces the multi-dimensional 
Burgers equation is widely used as a model of randomly 
driven Navier-Stokes equation without pressure [HI . [l2[ . 
[l3, 14, 15]. The three-dimensional form of equation has 
been used to model the formation of the large scale struc- 
ture of the Universe when pressure is negligible. Known 
as " adhesion" approximation this equation describes the 
nonlinear stage of gravitational instability arising from 
random initial perturbation [H, [H [H, [ll [13] . Other 
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problems leading to the multi-dimensional Burgers equa- 
tion, or variants of it, include surface growth under sput- 
ter deposition and flame front motion [U, [11] . In such 
instances, the potential tp corresponds to the shape of the 
front's surface, and the equation for the velocity potential 
tp is identical to the KPZ (Kardar, Parisi, Zhang) equa- 
tion d, [m, [H, [11]. For the deposition problem the ve- 
locity in multi-dimensional Burgers equation v = —'Vip 
is the gradient of the surface. The mean-square gradi- 
ent E{t) — ((VV'(x, t))^) = (v^(x, i)) is a measure of 
the roughness of the surface and may either decrease or 
increase with time. 

When the initial potential "00 (x) is a superposition 
of one dimensional potentials tpo,iixi), namely ipo{'x.) — 
J^i'^oA^i) and wo,i(x) = vo^i{xi), there are no in- 
teraction between the velocity component wo,i(x, i) = 
vo^i{xi, t) and evolution of each component is determined 
by one-dimensional Burgers equation. Before the de- 
scription evolution of fields in multi-dimensional Burgers 
equation we discuss now very short the evolution of basic 
types of perturbation in one-dimensional Burgers equa- 
tion [J, [1, [1, [1] , and compare the behaviour of "plane" 
orthogonal waves in 2-dimensional Burgers equation with 
different initial spatial scales. 

At infinite Reynolds number (y — > 0) the harmonic 
perturbation V{i{x) = fcoV'o sin fcox {iPq(x) — 4'ocoakox), 
is transformed at i ^ t^i = l/fcgV'o into saw-tooth 
wave with gradient dv/dx ~ 1/t and the same period 
Lq — 27r/fco. It's important that at this stage the am- 
plitude a = Lo/t and the energy E{t) = L^/12t^ doesn't 
depend on the initial amplitude. Thus if we compare 
the evolution of two components vo_i{xi,t) with equal 
potential ipo and different scales Li (Li << L2), the 
initial energy will be much higher for the component 
with smaller scale i?i(0)/£'2(0) = L\/ L\. But asymptot- 
ically we have inverse situation Ei{t)/E2{t) 
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For large but finite Reynolds number Reo = tpo/^v 
the shock fronts have a finite width ~ vt/L^ and at 
t 3> ^niReo we have a linear stage of evolution where 
v{x,t) = 4i/fco sin(fcoa;) exp(— i/fept). 

Continuous random initial fields are also transformed 
into sequences of regions (cells) with the same gradient 
dv/dx = but with random position of the shocks 
separating them. The merging of the shocks leads to 
an increase of the integral scale of turbulence Lit) and 
because of this the energy E{t) ~ L'^{t)/t^ of random 
field decreases more slowly than the energy of periodic 
signal. The type of the turbulence evolution is deter- 
mined by the behaviour of large scale part of the initial 
energy spectrum Eo{k) ~ a^fc". For n < 1 the initial 
potential is Brownian or fractional Brownian motion and 
scaling argument may be used [1, 0, H, [H, [H, H^. In 
this case the turbulence in self-similar and with integral 

2 

scale L{t) = (at) 3+1. For n > 1 the law of energy de- 
cay strongly depends on the statistical properties of the 
initial field H [2I [H, Hi, [s^, [Mi . 

For an initial Gaussian perturbation the integral scale 
L{t) and the energy of the turbulence E{t) 

E{t) = L'm/t^ (2) 

are determined only by two integral characteristic of 
the initial spectrum : the variance of initial potential 
oj,^ = hj^) and the velocity al = {vD 0, i, [32, [13, [13, 
[ssl . [3R [331 ■ Here Iq = jov is the integral scale of initial 
perturbation, and o^jl"^ = tni is the nonlinear time. Thus 
the energy of two components with equal initial potential 
variance and different scales Zo,ij ('0,1 '0,2) will have 
very large difference at t = 0; Ei{0)/'E2{0) = Ij/lj > 1 
and with logarithmic correction will be the same at large 
time Ei{t)/E2{t) ~ 1. For large, but finite Reynolds 
number Rcq = a^/2v the shock fronts have a finite 
width ^ vt/L{t) and due to the multiple merging of 
shocks the linear regime take place at very large times 
t ^ tniexp(Reo)/Reo. At linear stage the energy decays 
as Cf~^/^, where C ^ Lq exp(Reo)/Reo. 

The goal of the present paper is the investigation of 
the evolution of anisotropic regular structures and turbu- 
lence at large Reynolds number, when we have a multiple 
interaction of the spatial harmonics of the initial pertur- 
bation. We shown that we have local isotropization of 
the velocity and potential fields inside the cells. For the 
periodic wave we have the decay of frozen structure. The 
global structure of the random field is determined by the 
long correlation of initial field, and for the short cor- 
related field we have isotropization of turbulence. The 
other limit we consider in the paper is the old-age be- 
haviour of the field, when the processes of nonlinear self- 
action and harmonic interaction seems to be frozen, and 
the evolution of the field is determined only by the linear 
dissipation. 

The paper is organized as follows. In Section|TT]we for- 
mulate our problem and list some results about the so- 



lution of multidimensional Burgers equation in the limit 
of vanishing viscosity and its old age behavior. We also 
shown that we have local isotropization of the velocity 
and potential fields. In Section [Hi] we consider the inter- 
action of plane waves and evolution of periodic structures 
in 2-d Burgers equation. In Section HVl we consider the 
evolution of anisotropic multidimensional Burgers turbu- 
lence in inviscosid limit. We also discuss here the influ- 
ence of finite viscosity and long range correlation on the 
late stage evolution of Burgers' turbulence. 

II. MULTI-DIMENSIONAL BURGERS 
EQUATION, THE LIMIT OF VANISHING 
VISCOSITY AND THE OLD AGE BEHAVIOUR 

We shall be concerned with the initial value problem 
for the un-forced multi-dimensional Burgers equation ([1]) 
and consider only the potential solution of this equation, 
namely 

v(x,t) =~V^(x,t) . (3) 

The velocity potential '4>{'x.,t) satisfies the following non- 
linear equation 

^ = i(VV.)^ + .v2^. (4) 

The equation for the velocity potential tp is identical to 
the KPZ (Kardar, Parisi, Zhang) equation 0, [U, [11, 
which is usually written in the terms of the variable h = 
■ 4>. The parameter A has the dimension of length 
divided by time and is the local velocity of the surface 
growth. Henceforth ft,(x, t) has the dimension of length 
and is the measure of the surface's shapeness. In this case 
V = — Vi/' is the gradient of the surface. The roughness 
of the surface is measured by its mean-square gradient 

E{t) = ((V^(x, t)f) = (v2(x, t)) = ^'(^) ' (5) 

i 

Here the angular brackets denote ensemble averages or 
space averages for periodic structures. For the one- 
dimensional homogeneous field E{t) is the density of en- 
ergy and always decreases with time. At the initial stage 
of evolution and in the limit of vanishing viscosity the 
multi-dimensional Burgers equation is equal to the free 
motion of the particles. In Lagrangian representation the 
velocity of the particle V(t;y) is a constant. Here y is 
initial (Lagrangian) coordinate of the particle. In one- 
dimensional case the increasing of the length of elemen- 
tary Eulerian interval Ax = Ay -I- t/^V is compensated 
by the decreasing of length Ax = Ay — t/SV of the steep- 
ening interval and therefore the energy of the wave (the 
mean roughness of the curve) is conserved. After shock 
formation(s), the energy of the wave decreases with time. 
In multi-dimensional case the changing of elementary Eu- 
lerian volume depends also on the initial curvature of per- 
turbation and we don't generally have compensation of 
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steepening and stretching volumes. Thus for d > 1 the 
roughness of the surface, measured by its mean-square 
gradient E(t) (se e ([5])) may either decrease or increase 
with time [2 71 . |3l| . The increase of the mean-square gra- 
dient in the multi-dimensional Burgers equation (in con- 
trast with d = 1) is the result of this equation not having 
a conservation form. Nevertheless we will use the ex- 
pression "turbulence energy" for the value of E{t) and 
call Ei{t) = {(dtjj/dxi)^) — (vf) the energy of the i-th 
velocity component. 

Using the Hopf-Cole transformation [H, [s^ 



tp{x,t) = 2zylog[/(x,i) 



(6) 



one can reduce equation ([T]) to the linear diffusion equa- 
tion 



dU 
'dt 



vV'^U, [/(x,0) = J7o(x) =exp 



•00 (x) 

2iy 



(7) 



The goal of the present paper is the investigation of the 
evolution of regular structures and turbulence at large 
Reynolds number. For the not very large times we can 
use the the solutions of Burgers equation in the limit of 
vanishing viscosity. The other limit is the old-age be- 
haviour of the field, when the processes of nonlinear self- 
action and harmonic interaction seems to be frozen, end 
the evolution of the field is determined only by the linear 
dissipation. In this case we have the linearisation of Hopf 
Cole transformation 

In the limit of vanishing viscosity i' —^ use of 
Laplace's method leads to the following "maximum rep- 
resentation" for the potential velocity field [1, [l^ [ll] : 



'0(x, t) — max $(x, y, t) 
y 



$(x,y,t) ^ Vo(y) 



v(x,i) = 



y(x, t) 



(x-y)^ 

2t 



vo(y(x,i)) 



(8) 
(9) 



Here ?/'o(y) is the initial potential and vo(x) = — VV'o(x). 
In pop y(x, t) is the Lagrangian coordinate where the 
function $(x, y,t) achieves its global or absolute maxi- 
mum for a given coordinate x and time t. It is easy to see 
that y is the Lagrangian coordinate from which starts the 
fluid particle which will be at at the point x the moment 
t i. 

At large times the paraboloid peak in ([9]) defines a 
much smoother function than the initial potential "00 (y)- 
Consequently, the absolute maximum of <I>(x, y,t) coin- 
cides with one of the local maxima of ipQ (y) . In the neigh- 
borhood of local maximum yk we can represent the initial 
potential in the following form 



■00 (x) = 00, fc 1 



2L? 



(11) 



where Xi now are the principle axis of the local quadratic 
form describing the potential near the local maximum. 
At relatively large time the Eulerian velocity field v(x, t) 
in the whole space will be determined by the particles 
moving away from the small regions near the local max- 
imum of t/'o (x) : 



yj(x,t) = 



(1 + V^o,fct/i|) 



(12) 



Then, the Lagrangian coordinate y(x, t) becomes a dis- 
continuous function of x, constant within a cell, but 
jumping at the boundaries [1, [l^. The velocity field 
v(x, t) has discontinuities (shocks) and the potential field 
-0(x, t) has gradient discontinuities (cusps) at the cell 
boundaries. From (llOp . p^ it becomes clear that inside 
the cells the velocity and potential fields have a universal 
isotropic and self-similar structure : 



■0(x,t) = Vo(yfc) 



v(x,<) 



(x-yfc)^ 

2t 



Yk 



t 



(13) 



(14) 



One can see that due to the nonlinearity there's the lo- 
cal isotropisation of the velocity field in the neighbor- 
hood of the local maximum of "00 (x). The longitudinal 
component of the velocity vector v(x, <) consists of a se- 
quence of sawtooth pulses, just as in one dimension. The 
transverse components consist of sequences of rectangu- 
lar pulses. At large times the global structure and evolu- 
tion of the velocity and potential fields will be determined 
by the properties of local maxima 0o(yfe)- For the ran- 
dom field wall motion results in continuous change of cell 
shape with cells swallowing their neighbors and thereby 
inducing growth of the external scale L{t) of the Burgers 
turbulence. 

Let us now discuss the old-age limit of the solution 
of Burgers equation. Consider a group of perturbation 
(10) with the bounded initial potential (-00 (x)^) < oo assum- 
ing that 00 (x) is a periodic structure or homogeneous 
noise with rather fast decreasing probability distribution 
of the potential 0o. For such a perturbation in C/(x, t) 
we separate a constant component U : 



[/(x, t)^U + U{x, t) = U{l + u(x, t)) 



(15) 



Here u(x, t) = U{'x.,t)/U is the relative perturbation of 
field J7(x, t). Inserting into ([7|) we see that U docs 
not depend on time. Here C/o(x) and uo(x) are fields 
with zero mean value (on the period or statistically for 
noise). As times goes on, the viscous dissipation and 
oscillation (inhomogeneity) smoothing causes the ampli- 
tude (variance) of the field C/(x, t) to become less. At 
times when its relative amounts U is small in compare 
with U {\u\ ^ 1) the solution ^ is equal to 



v(x,i) -2uVUix,t)/U = -2j/Vu(x,t) 



(16) 
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As /7(x, t) and u(x, t) satisfy the linear diffusion equa- 
tion, then v(x, t) also at these times fulfills the linear 
equation. This testifies precisely to the fact that the evo- 
lution of the velocity field enters the linear stage. The 
accumulated nonlinear effects are described in this solu- 
tion by the nonlinear integral relation between the initial 
velocity field vo(x) and the fields [/(x, 0), U (|3l7p . and 
are characterized by the value | A-^ol/i' ~ Rcq. Here Ai/jq 
is the characteristic change in amplitude of i/jq, and Reg 
is the initial Reynolds number. 

From (fTB]) it is easy to get the well known result, 
that for Rcq 3> 1 the initial harmonic waves asymp- 
totically has also harmonics form, but with the ampli- 
tude not depending on the initial amplitude [U, At 
large initial Reynolds number the homogeneous Gaus- 
sian field vo{x) at the nonlinear stage transforms into 
series of sawtooth waves with strong non-Gaussian sta- 
tistical properties [1, [s^. Nevertheless at very large 
time, when the relation p6p is valid, the distribution of 
the random field v(x, with statistically homogeneous 
initial potential V'o(x) converges weakly to the distri- 
bution of the ho mog eneous Gaussian random field with 
zero mean value [27| . This stage of evolution is known 
as the Gaussian scenario in Burgers turbulence. In the 
absence of the long correlation of initial potential field 
the potential (and velocity consequently) have an uni- 
versal covariance function . But the amplitude of 
this function is nonlinearly related to the initial covari- 
ance function of the field "00 (x) and increases proportion- 
ally exp(i?eQ) with increasing of initial Reynolds number 
Rcq. When the initial potential has a long correlation ( 
((■0o(x)V'o(O)) = |x|-"F(x/a;) ,0 < a < 3 ) at large x we 
have conservation at linear stage as long correlation as 
the anisotropy of the field F{x/x) [13] • 

III. THE EVOLUTION OF PERIODIC 
STRUCTURES AND THE INTERACTION OF 
PLANE WAVES IN 2-D BURGERS EQUATION 

It was shown in previous section that at large time 
we have a cellular structure of the field with universal 
behaviour of potential and velocity inside each cell. The 
global structure of the field will be determined by the 
properties of local maximum of initial potential. 

Let us consider the evolution of periodic structure in 
2-dimensional Burgers equation 

ipoixi,X2) = 2V'o cos(fc2a;2) cos(fcia;i). (17) 

The evolution of this structure my be interpreted as the 
interaction of two plane waves with equal modulus of the 
wave number fco 

'>Po{xi,X2) — ipo cos(fciXi -I- k2X2) + tpo cos(fcia;i — ^22:2). 

(18) 

Here ki = koUi = 2tt/Ii, k2 = ^0^2 — 2-1: /I2 and n 
is unit vector with the component ni,n2- For the each 
of noninteracting plane waves the energy of i-th compo- 
nent is Ei^pit) = E{t)nf, where E{t) is the energy of 
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FIG. 1: Cellular structure of potential for ki S> ^2. 

plane wave and E{0) = For — > the plane 

wave is transformed at t ^ t^i — 1/fcQ^o into saw- 
tooth wave with gradient dv/dx — 1/t and the energy 
E{t) = TT^/B/coi^- For large but finite Reynolds num- 
ber we have a linear stage of evolution where v{x,t) = 
AiykQ sm(kQx) exp(— z^/cQi). 

Consider first the late stage of evolution of periodic 
structure of times t^tni- At this stage the velocity has 
the universal form in each cell (HH), where yk are the 
maximums of the initial potential (|17p . For the initial 
potential there are two sets of maximum of equal value 
corresponding the conditions cosfcixi — cosA;22^2 — 1 
and cosfcixi = cos/c2a;2 = —1. The shock lines (cell 
boundaries) of the velocity field are orthogonal to vector 
connecting the neighbor cell center and they are immobile 
and situated in the middle between the centers. Due to 
the symmetry of initial conditions the velocity field is 
symmetric over the point (^1/2,^2/2). Assume now that 
h < I2 and consider the velocity fields inside the region 
S: (xi G [QMI%X2 £ [0,^2/2]), see figure ??. 

The region S is divided by the shock line in the regions 
Si and 5*2 

s.,(o<.,<_ii.,,i!_La) . 

The center of the cell inside the region 5*1 is in the point 
(xi = 0,X2 = 0) and inside the region ^2 is in the point 
(xi = lx/2,X2 = h/2). Consequently for the velocity 
fields one can get 

Vi=Xi/t, V2^X2/t, xeSi, 

VI = {xi - h/2)/t, V2 = {X2 - l2/2)/t, X e ^2. (20) 

Thus on large time we have a frozen structure of the field 
with decreasing amplitude ~ For the energy of the 
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velocity component from expressions ([20)1 one can obtain (l22l) . For 'i^/'^k2 > fci > fc2 we have 



Ei{t) = 



E2{t) 



/2 

n 

12t2 
48*2 



1 - 



1 2 /cq n/'2i^ 



1 - 



1 + 



(21) 



From (|T9|) ■ ([20]) one can receive that for the very 
anisotropic fields (/i <C ^2), (j^i — 1,"2 ^ 1) the ve- 
locity component vi reproduce the behaviour of the ve- 
locity in one-dimensional Burgers equation, but the large 
scale component has now the period L ^ I2/2 instead of 
L — I2 for the initial perturbation. Let us compare the 
decay of the periodic structure (|17p . which is a super- 
position of two plane waves with the decay of the en- 
ergy component Ei^p of singular plane wave. For li <^ I2 
the energy of small scale decays as Ei ~ /ikQn\t^ ~ 
Ei,i = tt"^ n\ / ik'^t^ . The energy of large scale component 
E2 ~ 7rVl2fc^np2 > E^,2 = n'^nl/U^t'^ . 

Consider now the linear stage of the evolution of the 
periodic structure, when the Hopf-Cole solution is re- 
duced to the linear relation (|16p between the velocity 
field v(x, t) and the solution {/(x, t), u{x,t) of the linear 
diffusion equation ([7]). Using the relation exp {zcos6) = 
Io{z) + 2 J2m=i Irn{z) COS mO , whcrc Im(,z) are a modified 
Bessel functions, for the solution of this equation we have 
from 0,(1171) 



;7(x,i)= /2(Reo) 



00 

m— 1 



+2 



00 00 

EE 

n=l m— 1 



/o(Rco)/m(Reo) 

{cos[m(fcia;i -I- ^22^2)] 
-|-cos[m(fcia:i - ^2X2)]} 

/m(R.eo)/„(Reo) 

{cos[(n -I- m)kiXi -f (n — m)k2X2] 

^-~w{{n+m)^kl + (n-m)^kl))t 

+ cos[(n — m)kixi + {n + m)fc2a;2)] 



-v{{n~mykl + {n+mYk^)t 



} 



(22) 



where Reg = ipo/2i'. Here the first sum described the 
nonlinear evolution of two plane waves and the double 
summation - the interaction between the plane waves. 
From eq. (|22p we have that a constant component in eq. 
(Uni) is [7 = /o (Reo) and C7(x, t) = C/(x, t) ~ U. At hnear 
stage of evolution, when u(x, t) — U{x,t) /U ^ 1 we have 
from eq. (p)). p^ 



(23) 



tp{x, t) ^ ^ + Vj(x, i), -0(x, t) ~ 2i^-u(x, i). 



where ^ = 2v log U . The asymptotic behaviour of the 
potential (shape of the surface) V'(x) will be determined 
by the low index of the decaying exponent in solution 



Vi(x,i)~ 4i^(/i(Reo)//o(Reo)) 

X cos(fciXi)cos(fc2X2)e-''('=i+'=')* (24) 

and consequently the velocity component Vi decays like 
the j-th component of the velocity of plane wave 

?;i(x,t)= 4z.fci(/i(Reo)//o(Reo)) 

X sin(fcia;i)cos(fc2X2)e-''('=?+'=')* (25) 

W2(x,<)= 4!.fc2(/i(Reo)//o(Reo)) 

X sin(fc2a;2)cos(fciXi)e-''('=i+'=^)* (26) 

For the small Reynolds number this solution is equal 
to the linear solution of Burgers equation. But for 
3^/2^2 < ki the nonlinear interaction between the plane 
waves change the asymptotic evolution of the potential 
(shape of the surface). Now the leading term in eq. ((^ 
is along axis xi 

Vi(x,t)~ 2i/(/2(Reo)//o(Reo)) 

X cos(2fc2X2)e-''('*'=^')* (27) 

and is on the second harmonic along axis X2. Thus 



U2(x,i)= 4z.fc2(/2(Reo)//o(Reo)) 
X sin(2fc22;2) 



(28) 



and we have depressing of the modulation of potential 
along axis xi. The evolution of the gradient of surface 
along xi (velocity component vi is still determined by 
eq. (|25l) and they decay faster than the gradient over xi. 

Thus due to the nonlinearity which leads to the gen- 
eration of cross-wave numbers we have for the velocity 
component V2 at linear stage instead of initial spatial fre- 
quency ^2 the leading term at the second harmonic. This 
one is true even for the small Reynolds number. For the 
large Reynolds number Rcq we have /i(Reo)//o(Reo) — 1 
and the amplitude of vi,V2 do not depend on the ampli- 
tude of initial perturbation. 

Let us consider the transition processes of very 
anisotropic field when the angle between interacting 
plane waves is small ni ^ n2- Consider first a more 
general situation when plane periodic is modulated by 
large scale function Af (x) and the initial potential is rep- 
resented as 



-00^ (x) = Af (x)V-'o cos (fcixi 



(29) 



We assume that the function A/(x) characterized by the 
scales I/M,z and Lm,; ^ = l/27rfci. For the plane in- 
teracting waves p7)) Af(x) — 2cos(fc2a;2)- For the initial 
velocity field we have from eq. 



vi^q{xi,X2) ^ kiipo sm{kiXi)M{xi,X2), (30) 
t^2,o(a;i,a;2) = -■0oCOs(fcia;i)A4 (a;i,X2). (31) 
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In the limit of vanishing viscosity v the evolution 
of the velocity field is described by the equation (fTU]) 
and v(x, t) = vo(y(x, f)), where y(x, i) is Lagrangian 
coordinate from which starts the fluid particle which will 
be at at the point x the moment t. While LM,i ^ the 
velocity component V2 <C vi and at t <C tni.2 = L^n 2/V'o 
one can neglect the drift of the particles along X2 axis. 
In this case for the Lagrangian coordinates we get 



Xi{yi,y2,t) = yi +tkiipQsm{kiyi)M{yi,X2) , 

X2{yi,y2,t) = 2/2 • 



(32) 



Before t < t„i.i — 1/kfipo this solution v(x, t) = 
vo(y(x, <)) is single- valued. For t > tni,i we need to intro- 
duce the shock in multi-valued solution. In quasi-static 
approximation we assume that the evolution of the ve- 
locity component vi{xi, X2, t) is equal to the evolution of 
initial harmonic perturbation vi^q = Asin^kiXi) in one- 
dimensional Burgers equation. The amplitude of the per- 
turbation A — kiiljoM{xi^rn,,X2) (xi.m = hm) depends 
on the coordinate X2 as a parameter and is assuming to 
be a constant of a;i on the each period of the harmonic 
perturbation. 

Let's now consider the nonlinear stage of evolution 
^ni,i ^ i ^ ini,2 when the velocity component vi trans- 
forms into sawtooth waves. Consider the region when 
M{xi,X2) > 0. It's easy to see that at t:^ tni,i each pe- 
riod /i(m — 1/2) < xi < li{m+ 1/2) will be cover by the 
particles from small region near the point xi = lim and 
the solution of the equations xi = X{yi,y2,t), X2 — X2 
may be written as 



2/1 - xi> 



Xi — Xi ^ 



1 +tkl'lpoM{xi.m,X2) 



2/2 = X2- 



(33) 



The shocks are situated at the line xi^s — li{m + 1/2). 
From ([5^ 1551) one can receive for the velocity component 
at li{m - 1/2) < xi < li{m + 1/2) and for Af > 



Vi{xi,X2,t) 



V2{xi,X2,t) 



xi - yi,, 



t 



1 



cos 



1 -I- tkf^poM{xi.m,X2) 
-llJoM^^{xi,m,X2) 

kixi 



tkl^l;oM{xi^rn,X2) 



(34) 



From this equation one can see that for t 3> tni,i the 
velocity component vi {xi ,X2,t) is transformed into saw- 
tooth wave vi ~ (xi — yi^rn)/t like in one-dimensional 
case. It means that we have a fully depression of ini- 
tial amplitude modulation of this component. The ve- 
locity component V2 loss the periodic modulation over 
xi and V2.oixi, X2,t) ~ —ipoM^^{xi,X2) for positive M 
and V2fl{xi,X2,t) ~ il;oM^^{xi,X2) for negative M. The 
energy of this component increases twice in compare with 
the initial energy. 



FIG. 2; 
1^1(3:2,*). 



Evolution of 



For such wave the velocity field may be also represented 
in the form 

v(x, t) ^ V((x, t) + v,(x, t), Vi(x, t) = (v(x, t)) (35) 

where the brackets (...) means the averaging over period 
li . Here v; (x, t) is the large scale component and Vg (x, t) 
is the small scale component. We assume that the evolu- 
tion of small scale component Vc, may be described in the 
quasi-static approximation. The mean velocity v; has the 
scale in order of LM,i, and at stage t <^ tni,j\/ = L\.i ^/cTip 
one can neglect nonlinear distortion and dissipation of 
this component. Then from eq. ^ we have 

^^^ = -iv(v?(x,i))=-ivi?.(x,i) , (36) 

The integration of this equation over t give the evident 
expression for the coherent component 



v,(x,i) = -V(^(x,i)) 



(37) 



Here we have used the equation ^ for the potential 
ip{x,t). From eq. (|36|1 we have that the generation of 
large scale component is determined by the gradient of 
the energy of small scale component. For the periodic 
modulation Es{x,t) a.t t ^ tni,i does not depend on the 
initial amplitude. It means that at this times there are 
no generation of the large scale component. The gradient 
of mean potential (i/'(x,i)) = ■>Jjo\M{x)\ - ll/l&t at this 
time does not depend on t and from eq. (|37|) we get 



(38) 



vz(x,t) = -V(^(x,t)) = -^oV|Af(x)|. 



The amplitude of the small scale component is li/t, while 
the amplitude of the large scale component is in order 
tPo/Lma and it means that a,t t > LM^h/i^o the main 
part of energy is in the large scale component. The non- 
linear distortion of large scale component is significant at 
t > min{L\.j The future evolution of this compo- 

nent is strongly depend of the properties of modulation 
function Af (x). 

For the periodic structure ([T7)) when k\ S> ^2 we have 
that the plane wave i/jq cos(A:ia;i) is modulated by large 
scale function Af(x) — cos(A:ia;i) eq. (|29|) . The initial 
perturbation in this case is periodic structure with peri- 
ods li = 27r/fci , I2 = 27r/A:2 and li <^ l2- 
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FIG. 3: Evolution of 

V2ix2,t). 




FIG. 4: Evolution of the velocity components energies. 



Before the nonlinear distortion of large scale compo- 
nent t <^ tni,2 = h/f^ip the evolution of structure take 
place like in general case of modulated wave eq. ([29]) . 
The velocity component vi{xi, X2,t) is transformed into 
saw-tooth wave (Figure [21) and we have a fully depression 
of initial amplitude modulation of this component (Fig- 
ure [2]). The velocity component V2 loss the periodic mod- 
ulation over xi and V2,o{xi, X2,t) ~ ipo\cos{k2X2)\'x2- 
The period of this component is twice less than the initial 
period (Figure [3]) and the energy increases twice in com- 
pare with the initial energy. The behavior of the energy 
is shown on the Figured 



IV. EVOLUTION OF ANISOTROPIC 
MULTIDIMENSIONAL BURGERS 
TURBULENCE 

A. The intermediate stage of evolution of 
anisotropic random field 

In this section we will consider the intermediate stage 
of evolution of anisotropic random field in the two di- 
mensional Burgers equation. Let us assume that the ini- 
tial potential ^'0(2^1, X2) is random and strong anisotropic 



field with the spatial scales h ^ l2- The initial energy 
Ei{t) = {vf) — {{dip / dxi)'^) — u^/lf of the velocity com- 
ponent vi is in this case much larger than the energy 
of the large scale component V2 ■ We can introduce the 
nonlinear time of z-th component as ini,i = ^f/ci/;- For 
t ^ tl\\,2^ the drift of the Lagrangian particles in direction 
X2 is relatively small. Then one can assume 2/2 = X2 in 
eq. © and consider the one-dimensional problem with 
the initial potential ipo{yi^ X2,t), where X2 is a parameter. 

Due to the condition li <^ I2 the the first shock lines 
in the Lagrangian coordinates are on the points where 
dip/dyi has a minimum. In Eulerian space the are ori- 
ented primarily along the X2 axis end its length in this 
direction increase in time. For t ^ tni.i the velocity 
field vi{xi, X2,t) transforms to the sequence of triangu- 
lar pulses 



Vi{xi,X2,t) 



xi - yi,k{xi,X2,t) 

t 



(39) 



where yi^k{xi, X2-, t) are the coordinates of absolute max- 
imum of 19]) over yi with 2/2 = 2:2 ■ The shock positions 

xlj.{x2,t) are 



yi,fc+i + 2/1, fc 



+ Vkt 



(40) 



i^o{yiM{xi-,X2),X2) - 1po{yiM+l{xi,X2).,X2) 

yi,k+i{xi,x2) - yi.k{xi,X2) 



It means that at fixed X2 the interval x^ ^ < xi <, x^ i^j^-^ 
will be cover by the particles from small region near the 
Lagrangian point yi^k{xi,X2,t), and for the velocity com- 
ponent V2 we get 



V2{xi,X2,t) 



dlpQ{xi,X2). 



t:i=yi,k{^2;t) J 



dx2 

for xIj, <x< a:^+i . 



(41) 



The velocity V2{xi, X2,t) doesn't depend on xi between 
the shock-lines x^ f.{x2) and k+i^^'^)- "^^^ colhsion of 
the shocks in one-dimensional Burgers equation is now 
equal that at some point two adjacent shock lines a; J 
and x\ ^j^i touch each other. Then this point will be de- 
veloped into new shock lines x''^*^{x2) with the increasing 
in time length along axis X2 and which on its ends trans- 
forms into lines xl f^^xf f._^i. Thus at tn2.i ^ ini,i the 
velocity field has a cellular structure, the border x'l kix2) 
of the cells are describing by the equation (|40)) . the ve- 
locity component fi(x) has an universal structure (j39p . 
Velocity component V2 (x) inside the cell does not depend 
on Xl and along the axis X2 reproduced the behaviour of 
V2 along the line xi = 2/i,fc(2^2i t) cq- (|4ip . The evolution 
of the potential is plotted on Figure [5l 

The statistical problems of the velocity component 
vi at this stage are similar to the properties of one- 
dimensional Burgers turbulence. The integral scale Li{t) 
and the energy Ei (t) jt of the component v\ are 
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described by the expressions ([2]), where cr^ is the vari- 
ance of initial two-dimensional potential iJjq and Iq — li 
is the integral scale of vi component Li = a^ja^^. Due 
to the merging of the shock lines the integral scale of 
the turbulence along the axis x\ increases with time 
L\(t) ~ (t(Ti^)^/2 and at t ^ ^ni,2, when L\{t) ~ and 
E\(t) ~ i?2(0) we need to take into account the non- 
linear distortion along the axis X2- At i ^ tni,2 the 
potential and velocity fields have a universal isotropic 
and self-similar structure inside the cells: ea. fT^ . pi)) . 
The boundary of the cells on this stage degenerate into 
straight lines (planes, in three dimensional case). The 
multiply merging of the cell will leads to the establish- 
ment of statistical self-similarity and isotropization of the 
field. In the next section we will show how the statistical 
properties of the isotropic multidimensional Burgers tur- 
bulence are connected with the parameter of anisotropic 
initial perturbation. 

B. Isotropisation of the multidimensional Burgers 
turbulence 




10 20 30 40 50 60 



The statistical properties of the Burgers velocity field v 
(equation (fTO|) ) in the limit — > are determined by the 
statistical properties of the absolute maximums coordi- 
nate y(x, t) of the function ([5]). In one-dimensional case 
the problem of the absolute maximum is reduced to the 
problem of the crossing the random signal "i^^ix) by the 
non homogeneous function (x — y)^ jit + 11. The asymp- 
totic behaviour of the field at large t is determined by the 
maximum which amplitude is higher than the variance of 
the initial potential. That's why one can use some results 
of the theory of extremal processes [3, IH, liO] . In multi- 
demensional case the problem of the peaks statistic of the 
Gaussian field is rather well known for the isotropic and 
homogeneous field '4l'|. But for the Burgers turbulence 
we need to find the statistical properties of the absolute 
maximum of the scalar non homogeneous and anisotropic 
field $(x, y, i). In paper [i^l it was shown that at large 
t this problem is reduced to the problem of the finding 
of the statistical properties of extremum of random field 
?/'o(x) which value is much greater than their variance 

Let us assume the initial potential ^/'o(x) is random 
Gaussian field and it's correlation function may be used 
in the following form 

d 

(V'O(X)V'O(X + Z)) = Bv,(z) = G% n R^{Z^). (42) 



21^2 



(43) 



We assume also that the correlation function decreases 
rather fast at large distances; -B^,(|z| > 1st) — 0. Then 
the Gaussian initial field V'o(x) in the points |xi — X2I > 
1st is statistically independent. 




FIG. 5: Evolution of potential. (Initial potential and stage of 
isotropization) 



In the limit of vanishing viscosity we have "maximum 
representation" for the potential eq. In this solu- 

tion the velocity field v(x, t) eq. (fTO|) is determined by 
the coordinate y of the absolute maximum of function 
<I>(x,y,t). Let Q(-ff, AI4) denote the cumulative proba- 
bility and Wmax{H, AVjt) denote the probability density 
of the absolute maximum in elementary volume AT4 

QiH, AVk) = Probi<P <H,ye A14), (44) 
WmaAH, A14) = Q'^iH, AVk). (45) 
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Here we introduce the elementary volume A 14 which 
scale is much greater than l^t, but much smaller than 
the integral scale of turbulence L{t). The probability for 
the absolute maximum to be contained between Hi and 
Hi + AHi with the coordinate y(x, t) E AVk equals to 
that one for the absolute maximum Hi E {H, H + AH) 
to lie in AVk and to be less than H in magnitude for 
outer intervals AV^ 

Prob (y E AVk,H E [Hi, Hi + AH]) 

- W^aAH, Avk)AHQ{H,AVk). (46) 

Here we propose that the absolute maximums are statis- 
tically independent in the intervals AVk and AVk ■ The 
probability for coordinate y(x, t) to fall into AVk can be 
obtained by the integration of with respect of H 

Prob{y E AVk) = J W„,ax{H, AVk)Q{H,'EV k)dH. 

(47) 

After the integration of ([T7|) by parts we have 

Prob{Y E AVk) = J N{H, AVk)Q'H{H)dH, (48) 

where Q(H) is the integral distribution function of abso- 
lute maximum in the whole space. In Appendix in (4^ 
it was shown that for large H the integral distribution 
function 



QiH, AVk) = exp i-N{H, AVk)) 



(49) 



is determined by the mean number of extrcmum 
N{H, AVk) with value larger than H. 

Let us first consider the statistical properties of the 
extremum of the homogeneous random field ■00 (x). It's 
known that for the smooth fields the number of crossing 
some higher level asymptotically tends to the number of 
maximum and number of extremum. It means, that all 
the peaks above some high lever have only one extremum, 
which is the maximum of this peak. Thus we will consider 
first the properties of extremum of the field "00 (x). Using 
these properties of ^- function one can obtain for the mean 
number N{H; V) = {N^xp) of extremum with the value 
greater than H in some volume V 



N{H; V) = {^J^ diVi^o)\J{a,j)\E{i^o - H)d^J . (50) 

Here E(s) is a unit function, and J is the Jacobian of 
transformation 



J{a. 



dipoi'x.) 



(51) 



For the homogeneous field V'o(y) one can introduce the 
density of extrcmum as next{H; V) — N{H)/V. The den- 
sity of Uext in this case is determined by the joint proba- 
bility distribution function of the ^pQ, their gradient Vi = 



dipo/dxi and tensor = dtl'o/dxidxj. For the homoge- 
neous Gaussian field W^^^y^^a.^ = Wy{vi)W^„^aij{ipo, a^j) 
and from eq. (|50p we have for the density of extremum 

UextiH) = Wy{0) J dS J day J(ay)W^„,Q.^(s,ai,)- 

(52) 

For the Gaussian field the p.d.f. of the field ipo{x) and 
it's derivative are determined by the correlation func- 
tion of tpoipc) eq. (UU. In equation ([5^ we will in- 
tegrate over the conditional probability WconicLij / S) = 
Wia,j,S)/W^,{S) and will get 



n,xt^Wy{0) J dSW^„{S) J da,jJ{a,,)W,orMj/S). 

(53) 

Using the properties of Gaussian variables one can receive 
that the conditional expected value of is {en j) con = 
S{aijS) /a'^. In the problem of the Burgers turbulence 
at large time the asymptotic of riext of high value H is 
important. Thus in conditional averaging we have 



{J{a,j))c 



J{{aij)con) ^ W{ 



'■O.eff 



(54) 



Here we introduce the effective length lo.eff 



id 



(55) 



and take into account that (a|) = cr^/Zo^i- Finally we ob- 



tain from equations (|53p , (j54p for the density of extremum 
the following expression 



/■oo 

ncxt[H) = W„(0) / dSW^,{S)j^ 

Jh leff 

1 



{2TT) — l'^ff JH/a^ 



d-1 



H 



2d 



(56) 



Thus from equation (j56p one can receive that the mean 
number of the extremum of anisotropic field ip^ (x) is de- 
termined by some effective spatial scale ^e//, which is 
geometrical mean of spatial scale Zo,i- For the relatively 
large H the density of extremum (|56p is equal to the 
density of events that the random field 0o(x) is over the 
H. 

For the homogeneous field N{H;V) = Vnext{H), 
where next{H) is described by the equation ([5S|) . For 
the nonhomogeneous field, even in one-dimensional case, 
the expression for N is more complicated. We assume 
that the nonhomogeneous field is $(x) = ipo{'x.) — a(x) 
and a(x) is a smooth function in scale of ipo{'x.). Then in 
a quasistatic approximation one can receive for the mean 
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number of events N{H; V) that $(x) > ff in a volume 
V the following expression 



N{H;V) = / ne:,t{H + a{x))dV, 



(57) 



where next{H) is determined by the expression ()56|1 and is 
the density of extremum of the statistically homogeneous 
function V'olx;). 

At large time the paraboloid a — (x — y)^/2t in equa- 
tion ([5]) is a smooth function in the scale of the initial 
potential. Then for the mean number of maximums one 
can use quasi-static approximation eq. (|57p and 



Q{H) = cxp {-N^{H)) 
NooiH) = 



next {H+^^^^^]dV. 



(58) 
(59) 



Here NooiH) is the mean number of extremum of 
<i>(x, y,t) in the hole space with magnitude greater then 
H and next [H) is the density of the number of extremum 
of the initial homogeneous potential V'o(x) with value 
greater then H . For H the density next{H) is 

determined by the expression ([SS]) and we have 



1 



(27r)(d+i)/2id 



H 

cr,p 



eff 



hV'' 1 



ait 



d/2 



(.60) 



eff , 



In equation (j58[) we integrate over the infinite space, but 
the effective volume {<^'^t/ Hll^ ^)'^/'^ is determined by the 
paraboloid term in equation ( [55]) . Now the effective 
number of independent local maximum in initial pertur- 
bation is Nmax Ve//)''^^ ^^'^ incrcascs with time. 
When Njnax ^ 1 we can introduce the dimensionless po- 
tential h as follows 



H = ha^, h ^ ho{l + z/hl), 



(61) 



where ho = Hqg^ and Hq is the solution of the equation 
iV(ffo) = 1 

/ ^ ^ \ 1/2 

ho^d^/^ Mog ^^ (2^)i/rf ) ' (V'(x,i)) ^(T^/io- 

^" . . (^2) 
Thus we have a logarithmic growth of mean potential. 

The dimensionless potential has double exponential dis- 
tribution 

Q{z) = exp (- exp (-z)) , 
Qh{h) = exp(-exp(-(/i-/io)/io), (63) 

One can see that for t ^ t^i = l^^j/a^ we have 
Nmax ^ 1 and the integral distribution of absolute 



maximum is concentrated in narrow region AH/H ~ 
a^/H^ <C 1 near Hq- Using this fact one can get from 



(j48|) the probability distribution of the maximums coor- 
dinate 



where 



VF(y,x,i) 



1/2 



1 



exp - 



2L2(t) 



(64) 



-1/4 



1/d 



(65) 

is the integral scale of the turbulence. From equation 
(fT0| . ([64|l we see that the one-dimensional probability dis- 
tribution of the velocity field is Gaussian and isotropic. 
For the energy of each component we have 



-1/2 



(66) 

Thus for the anisotropic initial field there's the isotropi- 
sation of the turbulence. 

For the multi-dimensional Burgers turbulence the two- 
dimensional probability distribution, correlation function 
and energy spectrum where found in [3, [l6| using so 
called "cellular" model. In this model is assumed that 
in different elementary volumes the initial potential are 
independent and that the potential has a Gaussian distri- 
bution. In this model there's a free parameter A which 
is the size of the elementary cell. In the present work we 
consider a continuous initial random potential field with 
given correlation function (|42l) . ((43l) . The procedure for 
calculation of the two-point probability distribution func- 
tion is nevertheless similar to that one used in I, p. It's 
easy to show that for the two-point P.D.F. we have the 
same expression as obtained in ^] for the cell model, only 
the size of the cell A we used to change with the effective 
spatial scale leff- The effective spatial scale lo,eff ([55]) 
is determined by the scales l^^i of initial correlation func- 
tion (|42)) . For the two-point P.D.F. , correlation function 
and energy spectrum we also have the self-similarity and 
isotropisation at large times. In particular for the nor- 
malized longitudinal and transverse correlation function 
of the velocity field v — v/cr^_i we have 



Bll{x) = {vilV2l) 



dx 



Bnn{x) = - {vinV2n) = ^'(i), 



(67) 



(68) 



where x — x/L{t) and 

POO 



1 



2 ./-oo 5(i£±£,exp 



dz 



+ .g(%^)exp 
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ds. 



(70) 



It may be shown that the function P{x) is the probabihty 
of having no shock within an Eulerian interval of length 
xL{t). As far as potential isotropic field are concerned, 
the normahzed energy spectrum e(fc) is formulated via a 
one-dimensional spectrum eMN{k) of the transverse com- 
ponent. The energy spectrum E^{k,t) is isotropic and 
self-similar 



(71) 



At large wave number k the discontinuity initiation leads 
to the power asymptotic behaviour E{k,t) ~ k^^. At 
small wave number is also has the universal behaviour 



E(k,t) 



(72) 



which has do with the nonlinear generation of low- 
frequency component. In particular for the three- 
dimensional turbulence wc have E(k,t) ~ kH^/'^. For 
the large, but finite Reynolds numbers, the "shocks" have 
a finite width 5 ^ vt/L{t) and relative width increases 

/ \ 1/2 

slowly with time S/L{t) ^ l\og{a^t/l'^j:j)\ . Thus at 
very large time we have a linear stage of evolution. 



C. The linear stage of evolution of Burgers 
turbulence 

Let us now consider the linear stage of evolution of 
random field, when the potential ^/'(x, t) and the velocity 
field v(x, t) eg. pB)) are linearly connected with the solu- 
tion u(x, t) of the linear diffusion equation eq. ([7]). Here 
m(x, t) are the relative fluctuations of the fleld L/(x, t) eq. 
(fT5)l . Introduce the spectral density of the field u(x, t) as 



E^{k,t) = J B^{z,t)e<^^^dz, 



(73) 



where Buiz^t) — (u(z, t)w(0, i)) is a correlation function 
of relative fiuctuation field U (x, t) . The evolution of the 
spectral density £'„ (k, t) and variance cr^ (t) of u are de- 
scribed by the equations 



Eu(k,t) ^ Euo{k)e 



^„o(k) =i;„(k,0) (74) 



(75) 



For the homogeneous Gaussian initial potential ■(/'o (x) we 
have from eq. ([6]) 



"^-^^^^ ^ (^ / 



exp 



B^{z) 



e^C'^^dz , (76) 



where B^(z) is the correlation function of initial poten- 
tial ipQ (x) . The condition of Burgers turbulence entering 
the linear regime is t ^ tun, where tun is determined 
from the equation cr^(t;m) ~ 1 . From eq. (|74p one can 
see, that the old-age behaviour of the scalar field u and 
consequently the velocity field v will be determined by 
the behaviour of the energy spectrum i?„o(k) at small 
wave numbers k. When the correlation function of ini- 
tial potential ■00 (x) may be represented in the form ((42)) 
at small p and B^{\p\ > 1st) — 0, then from eq. ([75]) we 
have, that the spectrum Euo{k) at fc ^ Rbq /leff is fiat 
and 



Ko(k = 0) = 



1 



^^-p^(/e///Reo)^exp (Re^ 



(77) 

where l^ff is determined by equation (|55p . The spec- 
trum of the field u at large time is isotropic and has an 
universal form 



Eu(k,t)^Due 



-2uk''t 



(78) 



and consequently isotropic is the velocity field d, [T^, 
[27| . The energy of each component decays as Ei(t) ^ 

D„(l/t)-(''+2)/2. 

Consider now the case when the initial potential has 
a correlation function ([1^ at small x and has a long 
correlation at large x [27| 

B^(x) =4(|x|/?,„„3)-"Fb(x/|x|), a>0. (79) 

Here the function i^s(x/|x|) describes the anisotropy of 
correlation function at long distances. If a > d the evo- 
lution of the spectrum of u fea 175)1 and the velocity v 
will be the same as in absents of long correlation [23] . At 
< a < d the energy spectrum of initial potential has a 
singularities at small wave number 

£;^(k)=a^C„Jkr-'^^^B(k/|k|), 0<a<d, (80) 

where the function i^£;(k/|k|) is determined by the func- 
tion Fb(x/|x|) and describes the anisotropy of potential 
spectrum at small wave number. It was shown 27| that 
in this case we have the conservation of anisotropy at lin- 
ear stage and asymptotic behaviour of the spectrum of u 
is determined by the equation ([71)1 where 



Eu{k)^E^{\<i)/[^v) 



R6o^;"„„|k| 



a — d 



FE(k/k) . (81) 



For the velocity spectrum it means that it reproduced the 
initial spectrum of velocity at small wavenumber multi- 
plied by the exponential factor exp {—2iyk^t). These re- 
sults was formulated [27 1 fo r the correlation function of 
velocity fields. In paper 27| also was shown that asymp- 
toticly the velocity field has Gaussian distribution. But 
the transformation processes to the linear stage are not 
trivial and may be estimated on base of spectral repre- 
sentation. The behaviour of the spectral density of the 
field u at small wave number k is determined by the tail 
of correlation function B,i,{x.) ([75)1 and from eq. ([7S)l we 
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have that the spectrum is described by the equation ([8T|) . 
But with the increasing of the module of wave number k 
the power anisotropic spectrum transformed to the flat 
spectrum (1771) . The wave number kaf where this trans- 
formation take place may be estimated from the condi- 
tion that at fc ~ kaf these spectrum are the same order 
and we have 



lo'ng 



l/(d-a) 



(82) 



The condition of Burgers turbulence entering the linear 
regime is determined from the equation cr^(tzm) — 1- 
The main contribution in the variance cr2(i) eq. (j75p we 
have from the flat spectrum eq. (|77p . The wave number 

hin ^ {vtlin)~^/'^ is 



kun^{Reo/leff)e-''^"ld. 



(83) 



Thus for the ratio of this two critical wave numbers we 
have 



kaf /ki 



lin 



Re. 



(Q+2)/((i-Q) / hong 







l/(d-a) 



,-Ro^Q/d(d-Q) 

(84) 



and for the large initial Reynolds number kaf <C kun- 
From the equation cr^(t/i„) ~ 1 we have that tun — 
iniofi'Reo ^ exp (2ReQ) and is extremely large in compare 
with the effective nonlinear time inicff — ^eff/'^i'- 
easy to see from ([S^ *) that the time of "isotropisation" 
tiso — i^k^j is much greater then the nonlinear time 
tiso ~ tnicffe^^°o"/f'(''-") and this difference increaces 
when the index a is near the critical value a ~ d. 



and we have primarily the energy decay primarily of 
this component (see Sectior lTV Ap . After the time t is 
greater the of the nonlinear time of the the component 
with the largest scale ini ^ — max^gj/cr0 we have the 
isotropisation of turbulence in the scales in order of the 
integral scale of turbulence L{t) ([55)1 . In Section FlV Bl 
we consider the situation in absence of long scale cor- 
relation. But based on the results of one dimensional 
case [3^ we may suggest that there are no influence 
a long scale correlation on the evolution of the energy. 
Nevertheless we still have conservation of anisotropy at 
large scales |x| ^ L[t) ([7^ . In the spectral representa- 
tion we have the conservation of the inital spectrum and 
anysotropy at small wave number, but at fc ~ ks{t) ^ 
the initial spectrum trasformes into selfsimular spectrum 
([7T|) with the universal behaviour E{k,t) ^ k^'^^^'> eq. 
((72| at kL{t) < 1. Let us define an energy wavenum- 
ber kL{t) — L^^{t) ^ {ta.^)~^^'^, which is roughly the 
wavenumber around which most of the kinetic energy re- 
sides. Hence, the switching wavenumber ks{t) goes to 
zero much faster than the energy wavenumber. Taking 
into account the finite viscosity we have that after the 
i ^ tun — inioffRe,^"'^ exp (2ReQ) the nonlinear evolution 
of the spectrum is frozen and only the linear decays of the 
small scales is significant ([71)1 . The frozen spectrum of 
the velocity potential has a critical wavenumber kaf bel- 
low which the spectrum is anisotropic and reproduce the 
initial spectrum, and at A; > kaf is flat eq. (|77p . Thus 
at taf 3> t 3> tiin this part of the spectrum will play 
the dominate role in the evolution of the velocity spec- 
trum, consequently the energy of all velocity component 
are equal. At t 3> taf the spectrum of the velocity will 
be reproduced the small scale part of the initial velocity 
spectrum multiplied by the exponential factor (|74p . and 
we have finally the anisotropic field. 



DISCUSSION AND CONCLUSION 



Let us now discuss the evolution of the turbulence in 
presence of anisotropy at small ([42]) or large (|79p spa- 
tial scales. In initial perturbation the energy of the ve- 
locity component is cr^/^o,i and is greater for the small 
scales Zo,i- At the initial stage the scale of the turbu- 
lence in this direction increases faster then in others 
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